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We present results of lattice computations using overlap fermions on a twisted mass background. 
Nf = 2 full QCD gauge configurations have been produced by the ETM Collaboration with very 
light pions (down to less than 300 MeV), with small lattice spacing (a « 0.09 fm) and large 
volumes (V/a^ = 24^ x 48). By profiting of the good chiral properties of the overlap operator for 
the valence quarks, it is also possible to have a precise (and unquenched) determination of those 
physical quantities where the chiral properties are crucial. In order to have unquenched results, 
we match the valence quark mass with the sea quark mass. We also perform computations with 
different quark masses in order to simulate (partially quenched) Strange and Charm quarks. A 
typical application is the computation of Bk, for which we present first results. 
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1. Introduction 



Dynamical overlap simulations are extremely expensive. One interesting possibility is to use a 
different regularization for valence and for sea quarks. In fact, valence quarks are much less critical 
from the costs point of view (since they only appear in the final measurements, as in the quenched 
case), but much more critical from the point of view of the symmetries. In particular we can use the 
"twisted mass" (tm) regularization for the sea quarks, and the "overlap" (ov) regularization for the 
valence quarks. This so called Mixed Actions approach ^ is very promising, since it can strongly 
reduce (or even completely eliminate) the operator mixing problem and has the potentiality of 
delivering the most precise and cost effective results in the near future. 

Violations of unitarity by lattice artifacts, which are expected, can be studied analytically 
within ChPT. They may also take the form of {0{a^) suppressed) double poles, just like in Par- 
tially Quenched QCD, but a closer inspection suggests that these might be small in practice 
Moreover, since the exact (twisted mass) sea quark matrix is available, they can also be studied 
numerically. This is important in order to keep lattice artifacts under control. A first test is re- 
ported in [^]. Numerical simulations using a similar "mixed" approach has been reported by other 
collaborations also in this conference [Q, ^, ^, 0, ||, ^. 

In this proceedings we present our first physical results obtained with this approach. The 
present analysis, which is done with limited statistics, is mainly meant to check the set-up that we 
are using in order to decide on possible improvements. 

The outline of this work is as follows. In the next section we describe the detailed set-up of our 
computation. In section |[ we give physical results on the pion sector. In section ^ we discuss the 
computation of renormalization factors, which is done using the RI-MOM method and the Ward 
Identities. In section |5[ we comment on our preliminary computation of Bk- 



2. Details of the computation 

The gauge background that we use in the present work consists in the Twisted-Mass gauge 
configurations which have been produced by the ETM Collaboration [|l^, |ll]]. We summarize here 
the main features. We use twisted mass fermionic action at full twist with A^j = 2 degenerate quarks, 
tree level Symanzik improved gauge action at jS = 3.9 which corresponds to a lattice spacing 
a « 0.09 fm. The volume is V /a"^ = 24^ x 48. In the present study we consider a single value 
of the sea quark mass aj^ = 0.004 (the hghtest available), which corresponds to a pseudo-scalar 
mass niji « 300 MeV. As mentioned in the introduction, these first results are obtained with a low 
statistics of 54 independent gauge configurations. For more comments about the choice of this 
background for sea quarks we refer to JTO| , 11]. 



Valence quarks are described by the overlap operator [12|: 



D{m) = {p-—)D + m, 

D = -{\ + ^=), A = aDw-p, (2.1) 

where Dw is the Wilson Dirac operator and p is a parameter that we set equal to one, in order to 
optimize the locality properties of D [|3|]. Before applying the overlap operator we perform a single 
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HYP-smearing transformation [13]. The computation of the propagators is done with point-Uke 
sources chosen randomly on the whole lattice. The inversions are performed by computing exactly 
the lowest 40 eigenvalues and then using the SUMR algorithm [14, 15] with adaptive precision 
[jl6|]. Thanks to a multiple mass procedure [17], which can be extended to the SUMR solver [p^], 
we produced propagators for a wide range of bare masses down to am = 0.006 and covering the 
Strange and Charm range. This brought a negligible loss of precision at high masses. The cost of 
the computation of one full propagator is equivalent to the cost of producing a few independent 
gauge configurations. In order to understand whether the continuum limit is convenient in this 
approach, it will be important to check how the above cost ratio will scales when a ^ 0, at fixed 
physical volume. 

In a previous report [^] we discussed a wide range of tests performed on smaller lattices and 
we will not repeat them here. We only mention that the comparison of the scalar correlator shown 
in [||] was not repeated in the larger lattice, since the low-mode averaging 118]- that is necessary 
to have a clean scalar propagator - is rather expensive and we prefer to look at more physical 
quantities first. 



3. Results in the pion sector 

The first quantity that we consider is the pion mass, since this is also what we use to match 
the valence quark mass with the sea quark mass. This is shown in Fig. |T[ The horizontal Une 
(with tiny error-bars) marks the pion mass obtained in the "unitary" (tm-valence, tm-sea) set-up. 
From this comparison, the matching point is estimated to be (in the overlap bare quark mass) at 
am = 0.0075(10). The matching of one quantity implies of course that other quantities are only 
matched up to lattice artifacts. The hope of this approach is that these are not too large in physical 
quantities. 

The pion decay constant fj^ can be computed in a number of ways. The most interesting one 
is the one which does not rely on any renormaUzation factor: 

A = ^K0|P|7r)|. 

This can be compared directly with the tm-valence tm-sea result [|^, which is also 0{a) improved. 
In this approach f-^ turns out to be about 10-15 % larger than in [[To|], at the matching point, but 
also the error-bars are of the same order of magnitude, and therefore still compatible. It is clear. 



from this analysis, that some kind of noise reduction techniques as those employed in [10] would 
be important. 

It is also possible to compare our results for the pion masses and the pion decay constants 
with Chiral Perturbation Theory. The necessary Partially Quenched formulae have been computed 



in [|I9|] and the corresponding finite volume corrections in lpO|]. This comparison is shown in 
Fig. |2[ The dashed lines show the fit at finite volume, while the solid ones show the corresponding 
extr apolations at infinite volume. This gives a value of /o which is larger than [p^], as is clear from 
the considerations above. 
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Figure 1: Matching the valence-valence and the sea-sea pion masses. 
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Figure 2: Fit of the data against Chiral Perturbation Theory at finite volume (dashed Unes). The solid lines 
are the extrapolations of the curves at infinite volume. The pion mass is plotted in a way to make the presence 
of non linear corrections more evident. 



4. Renormalization constants 



The renormalization factors have been computed with the RI-MOM method []21[]. This is 
possible since the gauge configurations had been (Landau) gauge fixed before the computation of 
the propagators. 

It is important to note that the tree level overlap operator is different from the Wilson oper- 
ator and for p = 1 the difference is significant at high momenta (which are above the cutoff, but 
still important in the RI-MOM procedure). To take this into account we define the quark field 
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renormalization constant as: 



^ / N -1 ^AyvPvS Hp)] , . r ■ 2/ N / ■ 2/<^P\ n2 

Z^{^,g) = -i—^-Ll—^iJl (0{p)= sm\ap) + {sm\^)-p)^ 

48 w(p)i:^sin^ap^ L 2 



The other definitions are unchanged with respect to []21[]. We computed the renormalization factors 
for all bilinear fermionic operator and for some choices of four fermions operators. In general we 
find that the chiral extrapolation is very stable, although the plateaux are not always completely 
clear. As an example, we show in Fig. ^ (left) the plateaux for the renormalization factors of 
the Vector and Axial currents. In these cases we obtain, in the chiral limit, Zy = 0.98(5) and 
Za = 1.05(5), where the errors are only statistical. These can be compared with the renormalization 
factors obtained from the PCAC Word Identities. The relation between bare and PCAC quark 
masses is displayed on the right hand side of Fig. ^. From this Ward Identity one can derive 
Za = 1.19(3), where the errors are only statistical. 

The RI-MOM method can also be used to determine the renormalization factors of the four 
fermions operators. In particular, in the next section, we are going to use the renormalization factor 
of the operator i^a5=2 = [^Yv{^ — Ts)*^] [■^7v(l — Ys)^]- The RI renormalization factor can then be 
converted into the renormalization group invariant one using the anomalous dimension computed 
in [22]. This gives us Zg^^ = 1.48(3). The momentum dependence of Z|^(/i) and Zf^^ in the 
chiral limit are shown in Fig. ^ 



5. To- wards the computation of Bk- Comments and conclusions 

An obvious quantity which is particularly interesting in this approach is Bk, the Kaon bag 
parameter, which is related to the mixing of and by the expression: 

(^Va5=2(m)I^'') = yMifiB^(At) 

In fact a precise non perturbative determination of Bk would have a strong impact on the deter- 
mination of the associated CKM matrix elements. Moreover, it is only with an exactly chirally 
symmetric regularization that the operator ^a5=2 cannot mix with other operators (without need of 
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Figure 4: Left: The momentum dependence of Z^(/i) (bottom) and Z^^' (top) at the chiral limit. Right: 
The bare factor as a function of the pseudo-scalar mass. The Kaon mass corresponds here to am\ ~ 0.05. 



relying on any tuning procedure). Finally, we have now the possibility to remove the quenching 
errors. 

We computed Bk in a standard way employing the propagators described above. More pre- 
cisely we use the same procedure described in [p3[], although we use lighter quark masses. In par- 
ticular it was important to use the left hand current. Our results are shown in Fig. ^ and imply for 
the bare B-parameter B^^ = 0.66(7) and for the renormalization group invariant one Bk = 0.98(11) 
(errors are only statistical). Although the error-bars become very large at light masses, they are still 
reasonable at the Kaon mass, which is relevant for Bk- Nevertheless, some kind of noise reduction 



technique would be probably helpful and we are currently exploring those used in [ |1(]| ] 



Comparison with ChPT has been performed using the formulae in [T^, 20], and the inclusion 
of appropriate lattice artifacts can be done following the procedure in [^]. 
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